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ABSTRACT 

Only a handful of supernovae (SNe) have been studied in multi-wavelength from radio to X-rays, 
starting a few days after explosion. The early detection and classification of the nearby type lib 
SN 2011dh/PTF lleon in M51 provides a unique opportunity to conduct such observations. We 
present detailed data obtained at the youngest phase ever of a core-collapse supernova (days 3 to 
12 after explosion) in the radio, millimeter and X-rays; when combined with optical data, this allows 
us to explore the early evolution of the SN blast wave and its surroundings. Our analysis shows that 
the expanding supernova Shockwave does not exhibit equipartition (e c /eB ~ 1000), and is expanding 
into circumstellar material that is consistent with a density profile falling like i?~ 2 . Within modeling 
uncertainties we find an average velocity of the fast parts of the ejecta of 15,000±1800 km/s, contrary 
to previous analysis. This velocity places SN 201 ldh in an intermediate blast- wave regime between 
the previously defined compact and extended SN lib subtypes. Our results highlight the importance 
of early (~ 1 day) high-frequency observations of future events. Moreover, we show the importance 
of combined radio/X-ray observations for determining the microphysics ratio e /eB- 



1 Cahill Center for Astrophysics, California Institute of Tech- 
nology, Pasadena, CA, 91125, USA 

2 Department of Physics, Marquette University, PO Box 1881, 
Milwaukee, WI 53201. 

3 Homer L. Dodge Department of Physics & Astronomy, The 
University of Oklahoma, Norman, OK 73019. 

4 Astronomy and Astrophysics, Eberly College of Science, The 
Pennsylvania State University, University Park, PA 16802, USA 

5 National Radio Astronomy Observatory, P.O. Box 0, So- 
corro, NM 87801, USA 

6 Benoziyo Center for Astrophysics, Faculty of Physics, The 
Weizmann Institute of Science, Rehovot 76100, Israel 

7 Carnegie Institution for Science, 813 Santa Barbara St, 
Pasadena, CA, 91101, USA 

8 IPMU, University of Tokyo, Kashiwanoha 5-1-5, Kashiwa- 
shi, Chiba, Japan 

9 Department of Astronomy, University of California, Berke- 
ley CA 94720-3411, USA 

10 Computational Cosmology Center, Lawrence Berkeley Na- 
tional Laboratory, 1 Cyclotron Road, Berkeley, CA 94720, USA 

11 Dunlap Institute for Astronomy and Astrophysics, Univer- 
sity of Toronto, 50 St. George Street, Toronto M5S 3H4, Ontario, 
Canada 

12 School of Physics and Astronomy, Tel- Aviv University, Tel- 
Aviv 69978, Israel 

13 Department of Earth, Atmospheric, and Planetary Sciences, 
Massachusetts Institute of Technology, Cambridge, MA 02139, 
USA. 

14 Australian Astronomical Observatory, P.O. Box 915, North 
Ryde, NSW 1670, Australia 

15 Computational Physics Inc., 8001 Braddock Road, Suite 
210 Springfield, VA 22151-2110, USA 

16 Pontificia Universidad Catolica de Chile, Departamento de 
Astronomi'a y Astrofi'sica, Casilla 306, Santiago 22, Chile 

17 Space Science Institute, 4750 Walnut Street, Suite 205, 
Boulder, Colorado 80301 

18 Spitzer Science Center/Caltech, Mailcode 220-6, Pasadena, 
CA 91125, USA 

19 Astrophysics Science Division, NASA Goddard Space Flight 
Center, Greenbelt, MD 20771 

20 Department of Astronomy, University of Maryland, College 



Park, MD 20742 

21 Center for Research and Exploration in Space Science and 
Technology, NASA Goddard Space Flight Center, Greenbelt, 
MD 20771 

22 Space Telescope Science Institute, 3700 San Martin Drive, 
Baltimore, MD 21218, USA 

23 INAFCT, Osservatorio Astrofisico di Catania, Via S. Sofia 
78, 1-95123 Catania, Italy 

24 Supernova Ltd, OYV #131, Northsound Rd., Virgin Gorda, 
British Virgin Islands 

25 Department of Astronomy, University of Texas, Austin, TX 
78712, USA 

26 Naval Research Laboratory, Code 7213, Washington, DC 
20375-5320, USA 



2 



1. INTRODUCTION 

SN2011dh in the nearby galaxy Messier 51 was dis- 
covered on UT 2011 May 31.893 by A. Riou; detected on 
June 01.19 UT by the Palomar Transient Factory (PTF, 
Law et al. 2009, Rau et al. 2009); and rapidly spectro- 
scopically classified as a type lib supernova (Arcavi et 
al. 2011a; 2011b). 

The proximity of M51 (as in Arcavi et al. 2011b we as- 
sume a distance, d = 8.03±0.77 Mpc) motivated searches 
for the progenitor star. A putative massive star detected 
by the Hubble Space Telescope in pre-explosion images 
of M51 very close to the SN position was described by 
Maund et al. (2011) and Van Dyk et al. (2011). The 
interpretation of the progenitor is uncertain. Based on 
the rapid cooling of the expanding SN ejecta, Arcavi et 
al. (2011b) argue that the progenitor star must have a 
radius smaller than that of typical red supergiants or the 
supergiant progenitor of the well-studied Type lib SN 
1993J, see e.g., Weiler et al. (2007) (R< 10 13 cm). In 
contrast, the astrometrically coincident star is a F8 su- 
pergiant. Reconciling these constraints means that the 
supergiant is either unrelated to the SN or was a com- 
panion of the (now erstwhile) progenitor star. 

For core collapse supernovae the interaction of the 
blast wave with the circumstellar medium can generate 
detectable radio and X-ray emission. Additionally the 
fastest moving ejecta is related to the size of the progen- 
itor (the more compact, the higher the velocity; Cheva- 
lier & Soderberg 2010). These two properties motivate 
early radio and X-ray observations. Thus, we began a 
program of centimeter and millimeter wave observations 
of SN2011dh with EVLA and CARMA respectively. We 
initiated X-ray observations with the Swift Observatory. 
Following the end of our early monitoring a long-term 
program at EVLA was launched by A. Soderberg and 
collaborators. Their results are reported in a recent pa- 
per by Krauss et al. (2012). 

Soderberg et al. (2012) analyzed the Swift X-ray and 
two epochs (day 4 and day 17) of radio observations of 
SN2011dh, and found an expansion speed of ~ 0.1c and 
a mass-loss rate of 6 x 10~ 5 M Q yr -1 . With these in- 
ferences they concluded that SN2011dh is probably a 
type cllb SN, namely a SN that originated from a com- 
pact progenitor star as opposed to one that originated 
from an extended progenitor star (ellb); see Chevalier & 
Soderberg (2010) and references therein for explanation 
of these proposed sub-types. A similar result was found 
by Krauss et al. (2012). 

Here we present a full set of centimeter and millimeter 
wave observations (from day 3 to day 12 after explosion) 
as well as an analysis of the Swift and Chandra X-ray 
observations of SN 201 ldh. To our knowledge, these rep- 
resent the most comprehensive (essentially daily) set of 
pan radio (5 GHz to 100 GHz) observations obtained at 
early times of a core collapse supernova. As a result we 
are able to probe the circumstellar matter at smaller radii 
(wherein one can expect to find deviations from homol- 
ogous flows) and in principle directly infer the density of 
circumstellar matter. 

The paper is organized as follows. The millime- 
ter (CARMA) and centimeter (EVLA) observations are 
summarized in Ej2j In Sj3]we carry out a standard analysis 
for the radio observations assuming a synchrotron self- 



absorbed model. We then summarize the X-ray obser- 
vations (Sj4j and present a combined radio-|- X-ray anal- 
ysis (which incorporates inverse Compton scattering of 
optical photons to the X-ray band) in [J5] We present 
our conclusions as well as review the returns obtained 
from early millimeter wave observations of core collapse 
supernovae (^SJ). Such a review is timely given the immi- 
nent operation of the Atacama Large Millimeter Array 
(ALMA). 

2. CARMA AND EVLA OBSERVATIONS 

Starting June 4, 2011, observations were undertaken 
with CARMA, centered on either 107 GHz or 93 GHz 
with an 8 GHz bandwidth in the D array configuration. 
For flux calibration we used either Uranus (whose flux 
was bootstrapped from observations of 3C345) or Mars. 
The compact source J1153+495 was used as a phase 
calibrator. The CARMA data were reduced using the 
MIRIAD reduction softwar d 27 !. The log of observations 
can be found in Table Q] and a montage of CARMA de- 
tections is shown in Figure [TJ 

Observations with the EVLA were undertaken in the 
framework of our program "PTF Transients in the Lo- 
cal Universe" (PI: Kasliwal) as well as a ToO program 
to observe Type II Supernovae (PI: C. Stockdale). The 
array configuration was A. 3C286 served as the flux cali- 
brator. We used the AIPS software to reduce the EVLA 
data. The first EVLA detection of the SN was in both 
the K (22.46 GHz) and Q (43.46 GHz) bands at fluxes 
of 2.6 mJy and 6.95 mJy, respectively. The log of the 
observations can be found in Table [1] 

3. STANDARD SYNCHROTRON ANALYSIS 

The radio emission arises from relativistic particles, 
which are accelerated at the shock, gyrating in mag- 
netic fields, generated in the post-shock gas (Chevlaier 
1982; Chevalier 1998; Weiler et al. 2002). A large body 
of work confirms that (for most supernovae) the early 
radio emission can be described by a synchrotron self- 
absorption (SSA) model (Chevalier 1998). While internal 
free-free absorption might also be present at early times, 
the lack of a steeper spectrum at the optically thick part 
of the spectrum suggests that SSA is dominant (see also 
Appendix 1). Hence, throughout this paper we adopt 
Chevalier (1998) SSA model. 

Assuming a simple power law for the energy of the 
relativistic particles, the synchrotron emission for a self- 
absorbed source can be described by 

for frequency, v below v ai where v a is the frequency at 
which the optical depth from synchrotron self-absorption 
is unity, and 

S » K ^^A^ (p+1)/2 ^ (p - 1)/2 , (2) 

for frequencies above v a - Here, B is the the strength of 
the magnetic field, R is the radius of the blast-wave, D 
is the distance from the observer to the supernova, / is 
the volume fraction of the radio emitting region, and the 

27 http:/ /bima. astro. umd.edu/miriad/ 
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Fig. 1.— CARMA image cutouts of SN2011dh in the 3-mm band. The source at the center is SN2011dh. The image 
cutouts are 40" on a side centered around SN2011dh's most precise position, based on our EVLA observations, at Right 
Ascension=13:30:05.104 and Declination =47:10:10.915 (J2000.0). The positional error is +/-0.01 arcsec in each axis (see text). 



TABLE 1 

Summary of early radio observations of SN2011dh 



Day 


Frequency 


Flux 




[GHz] 


[mjy] 


4.21 


4.8 


< 0.17 


4.21 


7.4 


0.203 ± 0.036 ± 0.01 


4.25 


22.5 


2.6 ± 0.07 ± 0.26 


4.23 


43.2 


6.95 ± 0.17 ± 0.7 


4.08 


107 


4.55 ± 0.33 ± 0.46 


4.08 


230 


< 3 


5.01 


8.5 


0.455 ± 0.046 ± 0.023 


5.12 


22.5 


3.95 ± 0.07 ± 0.4 


5.10 


43.2 


6.47 ± 0.14 ± 0.65 


5.22 


107 


3.66 ± 0.35 ± 0.36 


6.92 


93 


2.52 ± 0.27 ± 0.25 


7.12 


8.5 


1.06 ± 0.03 ± 0.053 


7.20 


43.2 


6.42 ± 0.17 ± 0.64 


7.22 


22.5 


6.89 ± 0.06 ± 0.69 


8.91 


93 


1.84 ± 0.31 ± 0.18 


9.02 


22.5 


7.46 ± 0.04 ± 0.75 


9.02 


8.5 


1.58 ± 0.03 ± 0.079 


9.06 


5.0 


0.42 ± 0.03 ± 0.021 


9.12 


33.6 


7.49 ± 0.06 ± 0.75 


10.95 


93 


1.61 ± 0.30 ± 0.16 


11.01 


43.2 


3.19 ± 0.15 ± 0.32 


11.03 


22.5 


8.17 ± 0.05 ± 0.82 


11.98 


8.5 


3.15 ± 0.06 ± 0.158 


12.02 


5.0 


1.22 ± 0.03 ± 0.061 


12.91 


93 


1.06 ± 0.31 ± 0.106 


15.08 


93 


1.05 ± 0.17 ± 0.105 



Notes - Day is given in UT days in June 2011. The errors presented in the table represent the rms error from each image and a 
systematic calibration error (5% for v < 20 GHz and 10% for v > 20 GHz) that should be combined in quadrature. 

energy spectrum of the relativistic particles is given by a relevant case here; see §3.1): 
power law, N(E) = N E~p. 

The SSA model does not attempt to predict the abso- 
lute or even the relative fraction of the two key com- 15 -i/io / / \ -1 ^ 19 / S p \ 9 ^ 19 ( D \ 18 / 19 
poncnts, the energy in relativistic electrons and the i? p = 8.8xl0 / oB I — I I — J \y[pc) 
strength of the magnetic field. However, the minimum _ ' 

total energy will be achieved at equipartition, i.e., when f J cm (3) 

the energy in the relativistic electrons is equal to that \ 5 GHz / 

of the magnetic fields or / c b = 1 (see Readhead 1994). 

Here / b is the ratio of the energy in relativistic electrons an j 

to that of the magnetic fields. 

The measurement of the single-epoch SSA spectrum 
(both the optically thick and the optically thin parts) _ 4/lg / f \ - 4 / 19 /S p \~ 2/19 / D \ ~ 4/19 

can be inverted to yield the radius and the magnetic field B p = 0.58/ oB ( — — I ( — I ( — — I 
at that epoch (Chevalier 1998). For p = 3 (which is the V J V yJ \ivipcy 
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Here, S p is the peak fluaF^I. and v v is the peak frequency. 
We assume / = 0.5 (see Chevalier & Fransson 2006). 

Equations Q] and [2] are adequate to describe the broad- 
band spectrum at any given epoch, as long as the emis- 
sion is SSA. The dynamics of the supernova shell (which 
in turn depend on the velocity profile of the blast wave 
and the radial density distribution of the circumstellar 
medium) determine R and v a . These dependencies are 
generalized by allowing for power-law variations in key 
quantities (Weiler et al. 2002) and leading to the follow- 
ing equations: 

where the absorption expression describes an internal ab- 
sorption by material mixed with the emitting component 
and assumes planar geometry. The SSA optical depth is 
described by 

where both Kl and Kb are proportionality constants 
that can be determined by fitting the data, and S de- 
scribes the time dependence of the optical depth. 

3.1. Blast-wave physical parameters at single epochs 

We analyze observations from the first five epochs: 
2011 June 4, 5, 7, 9, and 11. During the first four epochs, 
data were obtained within a span of ~ 4 hours and thus 
are considered to be nearly simultaneous such that we 
assigned a central average time to all flux densities of 
a given epoch. We note that the last C-, and X-band 
observation were made on June 12 rather than on June 
11 and therefore we obtain the C-, and X-band fluxes 
for June 11 by interpolation. For each epoch we fit the 
EVLA-CARMA spectrum to the SSA model character- 
ized by p, S p and v v (Equations Q] and [2]) . We assume 
equipartition. For each epoch, the SSA model is well fit- 
ted by the data. The electron energy index, p, averaged 
over the five epochs we analyze is p = 3.00 ± 0.18. Next, 
substituting our fitted parameters into Equations |3] and 
|4]we derive the magnetic field (B) and radius (R) at each 
of our five epochs. 

As shown in Figure 2, the blast- wave radius, R oc t m 
where m = 1.14±0.24 while the magnetic field is propor- 
tional to i~ 1±0 - 12 . Theoretically, the value of m depends 
on the density structure of the blast wave. In particular, 
once the blast- wave decelerates, m = 1 — (1/??) where n is 
the velocity power-law index of the blastwave (R oc v~ v ). 
r\ can be described as a function of the progenitor star 
polytropic index, n, where r\ = (n + 1)/(0.19 x n). For 
red supergiants (convective stars) 1 < n < 3/2 leading to 
8.7 < i] < 10, while for blue supergiants and Wolf-Rayet 
(radiative) stars n — 3 leading to r\ = 7. Unfortunately, 
our data lack the precision to discriminate between these 
two possibilities. Since theoretically m < 1 we conclude 
that the mean velocity of the shock (computed as R/t 
where t is the time since explosion) is about 21,000 km 
s- 1 . 

28 The peak flux marks the transition of the spectrum from 
optically thick to optically thin. 
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Fig. 2. — The blast-wave radius (red circles) and magnetic 
field (red squares) estimates at 5 epochs. We assume an ex- 
plosion date of May 31.5, 2011. Relative to this explosion date 
the epochs are [3.7, 4.6, 6.6, 8.5, 10.5] days. The correspond- 
ing values are: R ps [0.64,0.77,1.19,1.46,2.01] X 10 15 cm and 
B Ri [5.9,4.8, 3.3, 2.7, 2. 1]G. Given the uncertainty in the peak 
fluxes, peak frequencies and p, we estimate the radius and mag- 
netic field uncertainties to be 11% and 9%, respectively. The 
curves are power law fits with R oc ^ 114 ±0- 24 (solid blue line) 
and B oc t~ 1± - 12 (dashed green line). 

Next we can place a lower limit on the total energy. In 
principle, in addition to the energy in the electrons and 
in the magnetic field some energy can be carried by pro- 
tons. We cannot estimate the proton energy. Thus the 
lowest estimate of the energy is given by the equipartition 
analysis and excluding that carried by protons (ions): 

£min = (1 + /cb) / (7) 

This yields (assuming / c b = 1 and / = 0.5) E m i n Rj 
[1.6, 1.8, 2.9, 3.8, 5.8] x 10 45 erg at epochs of June 4, 5, 7, 
9, and 11, respectively. 

3.2. Time- dependent solution 

We now fit the data to a comprehensive multi-epoch 
SSA model (as described by Equations [5] and [6]). The 
resulting fit parameters are a = —1.15, (3 = —0.96, Kl = 
453.43, Kb = 1.9772 x 10 5 , and d" = -3.42. This implies 
an electron power law distribution with p Rj 3 and R oc 
i - 94 which is consistent with what we found in the single 
epoch analysis we performed above. 

Assuming the CSM was formed by a constant wind 
mass-losfQ the CSM density structure will have the 
form pcsm oc r~ 2 . We next calculate the electron num- 
ber density which is given by n e = /eB^Efj 87r7 f meC ^ 
where 7, n is the minimum electron lorentz factor (Soder- 
berg et al. 2006). The resulting electron density is 

n e Ri 1.6 — 2.5 x 10 5 ( lo^cm ) 2cm_3 - I n order to es- 
timate the mass loss rate of the progenitor via wind 
prior to the explosion, M = 4irR 2 n e m p v w , an assumption 

29 We assume that the circumstellar matter is ionized. This is 
probably not an issue for blue super giants or Wolf-Rayet progen- 
itors. For red super giants, the strong UV flash should provide 
some amount of ionization (see Chevalier et al. 1982). 
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Fig. 3. — Radio emission as a function of time. In the up- 
per panel the result of the time-dependent fit (solid curves) is 
presented separately for each of the following frequencies u = 
[107,93,43.22,33.56,22.46,8.46,7.4,4.95] GHz from left to right, 
respectively. In the lower panel the time-dependent fit and its 
residuals for a frequency of v = 22.46 GHz is shown. The mea- 
surements presented in the lower panel were shifted to the latter 
frequency, assuming the fit parameters from §3.2. 



about the wind velocity, v w , has to be made. Following 
Chevalier & Fransson (2006), instead of assuming a wind 
velocity we will scale the mass loss rate by it and define a 
new parameter, A — M/4irv w . In the case of SN2011dh, 
4k5x 10 11 gem -1 , which is a factor of 3 smaller than 
the value derived by Soderberg et al. (2011) and Krauss 
et al. (2012). This value of A corresponds to a mass- 
loss rate of 10 -7 x (v w /10 km s _1 ) M Q yr -1 , consistent 
with a broad range of massive progenitor stars, from red 
supergiants to compact Wolf-Rayet stars. 

4. X-RAY OBSERVATIONS 

SN 201 ldh was observed with the Swift X-ray telescope 
(XRT; Burrows et al. 2003) in a series of observations 
beginning on 3.5 June 2011. We analyzed all observations 
taken through 23 June 2011 using the pipeline software 
of the UK Swift Science Data Centre at the University 



of LeicesteQ UK SSDC hereafter (Evans et al. 2007). 
To generate the X-ray light curve (see Figure [4]), we fix 
the source position at the known location of the SN and 
extract a light curve with half-day binning (43.2 ksec per 
time bin) using an aperture of 35.4". For background 
subtraction, we determine the time-average count rate 
in this identical aperture, 1.29 ± 0.14 cts ks _1 (0.3-10 
keV), from pre-SN2011dh Swift data. 

The SN X-ray emission appears to fade rapidly at 
early times, so we split the observations in two for pur- 
poses of X-ray spectral analysis, fitting data from the 
first three time bins separately from the rest. The data 
are then fit to X-ray spectral models in XSPEC v 12.7.0 
(Arnaud 1996), using a common Hydrogen column den- 
sity (Nh) across the full duration of the Swift obser- 
vations and allowing the power-law photon (r) index 
and time-average flux to vary between the two data-sets. 
The resulting fit yields the following spectral parameters, 
where uncertainties here and below are quoted at 90%- 
confidence: A H = l.l+H x 10 21 c m - 2 ; Ti = 1.1 ± 0.4; 



and T 2 = 1.9_ ; 4 . The results suggests a softening of the 
X-ray emission over the course of the Swift observations 
but note that the photon index overlaps for T — 1.5 at 
the 90%-confidence intervals. 

Prior to analyzing the resulting X-ray flux light-curve 
(Fig. we incorporate Chandra X-ray flux measure- 
ments from Pooley (2011) and Soderberg et al. (2011); 
these are converted to our chosen 0.3-10 keV band- 
pass using the authors' chosen spectral models and the 
WebPIMMS tooQ Analyzing the resulting X-ray light- 
curve, we find that it may be adequately characterized 
(in a x 2 sense) as a monotonic power-law decay emis- 
sion. The single power-law temporal decay index is 
ax = — 1.2 ± 0.2 over the interval of Swift and Chandra 
observations. The Chandra data confirm the absence 
of significant contaminating source emission within the 
Swift XRT aperture (Pooley et al. 2011; Soderberg et al. 
2011). 



5. COMBINED RADIO-X-RAY-OPTICAL ANALYSIS 

The X-ray photon index is 1.1 (early times) and 1.9 
(late times). The high frequency radio spectral index 
is —1 which corresponds to a photon index of 2. The 
X-ray emission is elevated by a factor of 50 relative to 
extrapolation of the SSA spectrum. Thus we conclude 
that the X-ray emission is not a result of synchrotron 
emission. Moreover, as already shown by Soderberg et 
al. (2012), thermal free-free X-ray emission is ruled out. 
Another explanation for the X-ray emission is that it 
arises from inverse Compton (IC) scattering of optical 
photons H 

Assuming that the X-ray emission arises from IC, the 
SSA and IC formulation of Katz (2011) can be used to 
infer R and B. For p = 3 



RozD(r slc F t ) 1/10 (S„, abs f 5 (!/)- 



(8) 



30 UK Swift Science Data Centre: http://www.swift.ac.uk/ 

31 WebPIMMS tool: http://heasarc.gsfc.nasa.gov/Tools/w3pinuns.html 

32 The SSA model yields B of a few Gauss. A v p of 30 GHz 
would require a Lorentz factor, 7 ~ 33. These electrons can inverse 
scatter SN optical photons to the X-ray band. 
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Fig. 4.— X-ray light-curve of SN2011dh (0.3-10 keV) as 
determined from Swift XRT observations (black diamonds) 
and the high-level analysis software at the UK Swift Science 
Data Centre (Evans et al. 2007). Error bars are la. Data 
have been binned to half-day intervals, with bins having less 
than 5% coverage dropped from analysis, and converted from 
counts to flux using the results of our X-ray spectral fits. 
The two red squares are Chandra measurements converted to 
(0.3-10 keV) band. A fit to the temporal decay is indicated 
as a single power-law decay (dotted green line). Counts-to- 
flux conversion factors (Swift observations) for the two spec- 
tral epochs (the first 3 bins and the rest of the bins) are 
ci, = 8.3 x 10" 11 erg cm" 2 s" 1 ct" 1 and c 2 , a = 5.3 x 10" 11 
ergcm~ 2 s" 1 ct" 1 , respectively, with respect to 0.3-10 keV 
absorbed flux, and ci = 8.8 x 10" 11 erg cm -2 s" 1 ct" 1 and 
C2 = 6.7 x 10" 11 ergcm~ 2 s _1 ct -1 , respectively, for 0.3- 
lOkeV unabsorbed flux. The latter are the counts-to-flux 
conversion factors we use to generate the X-ray flux light- 
curve for SN 201 ldh. 



B = 0.48 



r sic F t 



10 _12 erg cm _2 sec _1 



2/5 



mjy 



-2/5 



G, 



(9) 



.10 GHz, 

where r SIC = vS VtSyn / '(vS v jc), S„j C is the IC flux, 
Su,syn is the synchrotron flux, S u ^ a bs is the SSA flux, 
and F t is the optical flux. 

Applying the Katz (2011) equations to the optical (Ar- 
cavi et al. 2011b), X-ray and radio measurements on 
June 7, 9 and 11 yields B » [0.77,0.68,0.48] G. The es- 
timated error (obtained from application of the rule of 
error propagations) is m 15%. These values of B are 
w 4.2 smaller than the corresponding estimates obtained 
from an equipartition analysis ( £|3.1[) . Nominally, by us- 
ing equation @] again, / c b ~ 500 to 1700 with 10 3 as 
a reasonable mean (which we adopt). However, as can 
be seen from Equation [3] this only results in a 30% de- 
crease in the value of R, relative to that obtained from 
equipartition analysis (333J- Thus, the mean velocity 
is R/t 1.5 x 10 9 cms _1 . The electron density is now 
higher by a factor of ~ 50 and therefore the mass-loss 
rate is higher by a factor of ~ 20 (for a fixed wind veloc- 

ity)- 

6. CONCLUSIONS & RAMIFICATIONS 

In this paper we present the earliest millimeter- and 
centimeter-wave monitoring observations of the Type 
lib supernova SN2011dh in the galaxy M51, between 
four days and 12 days after the event. Using the wide 



frequency coverage of the EVLA (4 to 43 GHz) and 
CARMA (100 GHz) we measure the key parameters of 
the synchrotron self-absorbed spectrum. The radio ob- 
servations were accompanied by extensive X-ray observa- 
tions by Swift. The X-ray emission, argued to be inverse 
Compton scattering of the SN optical photons by the 
relativistic electrons that produced the radio emission, 
combined with the radio observations, allow us to relax 
the equipartition assumption and track the radius and 
circumstellar density. 

We infer a mean velocity of R/t rs 1.5 x 10 9 cms" 1 , 
for the SN Shockwave, with an uncertainty entirely dom- 
inated by the limitations of the theoretical framework. 
Within the framework of the model the velocity is 1.32 
to 1.68 xl0 9 cms" 1 . This value range is larger than the 
~ 10 9 cms _1 expected for an extended progenitor (red 
supergiant) but smaller than the ~ 3 x 10 9 cms _1 ex- 
pected for a compact progenitor. This may have impor- 
tant implications for the evolution leading to the forma- 
tion of the progenitors of SNe lib: if SNe lib are split into 
two distinct, well-separated classes, that would suggest 
they may arise from two different evolutionary scenarios. 
On the other hand, if additional intermediate objects are 
found , that may suggest a continuum of objects between 
ellb and cllb, a single progenitor class may be favored. 

We use the joint SSA+IC model analysis to infer 
the evolution of radius with time, R cx t m . We find 
rn = 1.14 ±0.24. The error bars are too large to directly 
constrain the physical nature of the envelope of the ex- 
ploding star (radiative versus convective). Assuming no 
additional systematics the combination of the early mea- 
surements with observations over several months may 
provide the means to determine the deceleration phase 
of the shock-wave in a more precise way and indirectly 
may shed more light on the density profile of the outer 
layers of the progenitor. 

Soderberg et al. (2012) infer a much higher shock- 
wave velocity, 3 x 10 cms , when assuming equipar- 
tition. As noted in §1 their analysis primarily rested 
on the X-ray data and optical dat43 (since the radio 
data was quite sparse, day 4 and day 17). The lim- 
itations naturally propagate to the robustness of their 
inferences. This caveat notwithstanding, these authors 
find f e B = 30, a factor of ~ 30 lower than the value we 
find, suggesting a Shockwave velocity of 2.5 x 10 9 cms" 1 . 
At least within the framework of the SSA+IC model we 
believe that our inference is more robust, thanks to our 
comprehensive radio and X-ray data sets. Moreover, as 
mentioned above, there is a factor of 3 difference between 
the derived equipartition mass-loss rate value in our anal- 
ysis compared to Soderberg et al. (2012) The origin of 
this difference is in the way this value is calculated and 
in the assumptions made. While Soderberg et al. (2012) 
assume e e = 0.1 wc do not make such an assumption. 
Our measurements suggest e ~ 0.3 which accounts for 
the different mass-loss rate values. 

Another set of comprehensive radio data was presented 
by Krauss et al. (2012) for days > 20 after explosion. In 
their analysis they assumed equipartition only and found 
an average velocity of « 25000 km/s. Furthermore, as- 

33 The analysis uses the method of Chevalier &; Fransson (2006) 
which describes the IC emission as a function of optical luminosity, 
f e B, shock- wave velocity, mass- loss rate, wind velocity, and time. 
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suming e e = 0.1, they found a mass-loss rate 3.5 times 
greater than what we find. In view of these differences 
we analysed the data presented in Krauss et al. (2012) in 
the same manner that we analysed the early data set we 
present above. Our analysis suggests an average shock- 
wave velocity of 18, 500 ± 2000 km/s at later times. Our 
derived lower velocity can be explained by a few fac- 
tors. First, our fitting method allows the electron energy 
power-law index, p, to vary. We find an average value of 
p = 2.8, while Krauss et al. (2012) keep the power-law 
index constant at p — 3. This leads to slightly differ- 
ent values of the peak flux and frequency which in turn 
leads to a lower value of the velocity. Another factor 
that contributes to the difference between Krauss et al. 
(2012) and our velocity value is the different coefficient 
used in Equation 3. While we are using the coefficient 
from Chevalier et al. (1998), Krauss et al. (2012) use a 
coefficient which is larger by 20%. When using Krauss 
et al. (2012) values for the peak frequencies and fluxes 
in Equation 3 using our coefficient, we find an average 
Shockwave velocity of 21,000 km/s which is consistent 
with our results given our errors. 

In light of the above uncertainty in the radius coeffi- 
cient and in the fitting method (constant vs. varying 
electron power-law index p), we note that the uncer- 
tainty in derived values of the Shockwave properties is 
greater than the error measurements alone. The differ- 
ent uncertainties, combined, can be as high as 40 — 50%. 
Therefore, any conclusion based on the absolute derived 
value of the Shockwave properties such as its velocity, is 
weakened. The time evolution of these properties is less 
sensitive to the above uncertainties and therefore may 
provide a more robust diagnostic. 

On another matter, thanks to broad banding, improved 
receivers, and a flexible correlator, the Expanded VLA 
(EVLA) has sensitivity gains ranging from factor of 2 to 
10 in the 1-40 GHz band. Separately, there have been re- 
lentless and steady improvements in the continuum sen- 
sitivity of millimeter wave arrays (CARMA, PdBI) with 
the ALMA now offering an order of magnitude increase 
in sensitivity in the millimeter and sub-millimeter bands. 

These great improvements offer powerful diagnostics 
in two ways. First, both synchrotron self-absorption and 
free-free absorption (from the circumstellar medium) de- 
pends strongly on frequency. High frequency observa- 
tions can be sensibly undertaken at very early times. 
(EVLA bands are usually self-absorbed at, say, day 1). 
Thus, early millimeter wave observations can probe the 
fastest moving ejecta (since in a homologous flow the 
fastest moving ejecta is at the greatest radius). Next, 
the intrinsic SSA spectrum is modified by the external 
free-free optical absorption. This is best measured by 
comparing lower frequency measurements to higher fre- 



quency measurements. A clear detection of free-free ab- 
sorption gives a model independent measure of the den- 
sity of the circumstellar matter. 

These two diagnostics motivated our CARMA+EVLA 
effort. The toy model given in Appendix A shows that a 
future SN, such as SN2011dh, if observed at, say, day 1, 
would allow us to meet at least one of these two goals. 
Given that surveys such as PTF are now moving to even 
faster cadence we would expect these two goals to be 
realized in the very near future. 
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APPENDIX 

WHY EARLY RADIO (CM AND MM) OBSERVATIONS OF SNE ARE IMPORTANT 
The free-free optical depth is 

T ff = 3.3 x IO-^-^i/q^EM (Al) 

where T = 10 4 T4 is the electron temperature (in degrees Kelvinj^3, i/qhz is the frequency in GHz and EM is the 
emission measure, the integral of r? e along the line of sight, and in units of cm -6 pc. 

34 We use the convention of X n = X/10 n where it is assumed, 
unless explicitly specified, that the units are CCS. 
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For a star which has been losing matter at a constant rate, M, the circumstellar density has the following radial 
profile: 

M 

p{r) = n(r) » = i^: (A2) 

where v w is the wind velocity, r is the radial distance from the star, n is the particle density and fj, is the mean atomic 
weight of the circumstellar matter. Thus n(r) oc r~ 2 . We assume that the circumstellar matter is ionized. This is 
probably not an issue for blue super giants or Wolf-Rayet progenitors. The strong UV flash for red super giants provide 
some amount of ionization and so a specific check needs to be done for such progenitors. The emission measure from 
a radius, say, r% to infinity is 



EM = 



(9^=in^(ri)- 3 (A3) 



where n* is the density of electrons at radius r„. We choose the following normalization for these two quantities: 
n = 10 6 ng and r = 10 15 ri5. The corresponding expression for the emission measure is 



EM = 1.08 x 10 8 n?« r 



.is (— ) 3 (A4) 



* Mis- 
using the above equation for EM and substituting r 15 = 0.0864 x vgt^ the free-free optical depth is 

r ff = 5.5 x 10 4 nl 6 rt 15 v^ 3 tfT^- 35 ^ (A5) 

We use SN2011dh as an example and set vg = 1.5, r»i5 = 1, and n»6 = 0.2. In this case the free-free optical depth 
is Tff = 6.5 x 10 2 T 4 1 ' 35 ^q^ z 1 . Setting the LHS to unity and assuming T4 = 1 yield the run of the frequency at which 
qualitatively the optical depth is unity: 

^GHz(ff) = 21.8*7 1 - 43 (A6) 

The SSA optical depth is described by 

where both Kl and Kb are proportionality constants that can be determined by fitting the data, and 8 describes 
the time dependence of the optical depth. In §3.2 we found the following parameters for SN2011dh: a — —1.15, if 5 = 
1.98 x 10 5 , and 5" = -3.42. The SSA optical depth in this case is 

r SSA = 7 x IO^g&V- 43 (AS) 

Setting tssa to unity yields the run of the peak SSA frequency as a function of time 

^ghz(SSA) = Ultf 937 (A9) 

Armed with equations A5 and A8 with the fitted parameters we found in §3.2, we can now look into the value of 
even earlier mm-observations (see Figure 5). If had we observed SN2011dh on day 1 the EVLA observations would 
have shown clear signature for free-free absorption. We also note that rg is proportional to n*r* whereas the SSA 
optical depth is a different function of and r*. Thus we can, in principle, obtain a different measure of and r*. 
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Fig. 5. — Theoretical radio spectra of SNe including synchrotron-self absorption and with (solid lines) or without (dashed lines) free-free 
absorption. Blue and red lines are for days 1 and 4 after explosion, respectively. The above spectra were produced using Equations 5, A5 
and A8 with the fitted parameters found in §3.2 for SN2011dh. 



